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Abstract 

In this paper a geometric method based on Grassmann manifolds and matrix Ric- 
cati equations to make hermitian matrices diagonal is presented. We call it Riccati 
Diagonalization. 

1 Introduction 

In this paper we consider a finite dimensional quantum model, so its Hamiltonian is a (finite- 
dimensional) hermitian matrix. In order to solve the model we want to make the Hamiltonian 
diagonal. 

Although we have a standard method for the purpose, to perform it explicitly is very 
hard (maybe, almost impossible). Let us make a brief introduction, see pQ, [2J. 
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Let if be a hermitian matrix, namely 



H e{H e Af(n; C) | if* = H}. 
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When we want to make H diagonal Elementary Linear Algebra shows the following diago- 
nalization procedure [A] =>■ [B] =>■ [C] : 
[A] First we calculate eigenvalues of H, 



= \XE — H\ — X n — tiHX n - L + ■■■ + (-l) n det H 



(2) 



There are n real solutions although it is almost impossible to look for exact ones, so let those 
be {Ai, A 2 , ■ • • , A n }. 

[B] Next we find eigenvectors {|A)} corresponding to eigenvalues 



H\Xj) = Xj\Xj) and (Xi\Xj) = S{ 

for all 1 < i, j < n. It is also almost impossible to carry out. 
[C] Last by setting 

f/=(|A 1 ),|A 2 ),---,|A n )) 



(3) 



we finally obtain 



H = UD H tf 



(4) 



(5) 



where D H = diag(Ai, A 2 , ■ ■ • , A n ) is the diagonal matrix. 

The procedure is standard, while to carry out it completely is another probleno Even 
if n = 3 it is very hard. In fact, for the hermitian matrix 



/ 



H 



hi a (3 



\ 



G H(3;C) 



a h 2 7 
V P 7 h j 

carry out the diagonalization. As far as we know it has not been given in any textbook on 
linear algebra. 



1 in Japan it is called a pie in the sky 
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Note that the characteristic equation /(A) = \XE — H\ is given by 

/(A) = \3-(h 1 + h 2 + h 2 )\ 2 + (h l h 2 + h l h 3 + h 2 h 3 -\a\ 2 -\f3\ 2 -\ 1 \ 2 )\ + 
\l\ 2 hi + \/3\ 2 h 2 + \a\ 2 h 3 - hih 2 h 3 - a/3>y - a/3^. 

To look for exact solutions by use of Cardano formula is not easy (for example, try it by use 
of MATHEMATICA or MAPLE). 

Therefore we present another diagonalization method based on Grassmann manifolds 
and matrix Riccati equations. 

2 Riccati Diagonalization 

In this section we show a new diagonalization method. Before it let us explain the idea with 
simple example for beginners. The target is 

(hi a \ 
)eH(2;C). (6) 
a h 2 J 

To diagonalize H above we consider a matrix 

i / 1 - A 

U = U(z) = , for ;£C. (7) 

\z 1 j 

It is easy to check U ] U = UU ] = 1 2 and \U\ = 1, so U is special unitary (U G SU(2)). 
Namely, U is a map 

[/ : C(C5 2 ) ^ SU{2). 

This map is well-known in Mathematics or Mathematical Physics. 
The calculation WHU gives 



U*HU 



1 z \ I hi a \ I 1 —z 
1 ^ \ z ? \ -z 1 / \ a ho I \ z 1 



hi + az + az + h 2 \z\ a — {hi — h 2 )z — az 2 
1 + \ z \ V a — (hi — h 2 )z — az 2 h 2 — az — az + h\\z\ 



so if we assume the equation 

a — (hi — h 2 )z — az 2 = -<=>■ az 2 + (hi — h 2 )z — a = (9) 
we have the diagonal matrix 

hi+otz+az+h2\z\ 2 



U^HU = 1+121 a . (10) 

hi — az— otz+h\\z\ 



From this two eigenvalues are obtained by 



l 1 2 I 1 2 

/ii + az + az + h 2 \z\ h 2 — az — az + hi\z\ , 

1 = 1 , I |2 ' A2 = 1 , I [2 ^ 1:L ) 

1 + Z 1 + Z 



under the equation ([2]), whose solutions are easily given by 



-(hi-h 2 )±^(h 1 -h 2 f + A 



a 



2 



2a 



(12) 



A comment is in order. The equation (Q is a special version of (generalized) Riccati 
equations. 

As shown in the example our diagonalization method is different from usual one (a kind 
of reverse procedure). Let us state our procedure. 

Riccati Diagonalization 

[A] For H G H (n; C) we prepare a unitary matrix U = U (Z) G U (n) where Z is a parameter 
matrix and calculate 

U ] HU EE W = (Wij). 

[B] We set 

Wij = for 1 < j < i < n 
and solve these simultaneous equations (a system of Riccati equations) to determine Z = 

(Zkl)- 

[C] We finally obtain the diagonal matrix 

W = diag(wn, w 2 2, ■ ■ , w nn ) 
where each component is an eigenvalue of H under [B]. 
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3 General Case 

In this section we consider a generalization of the example in the preceding section. See [3], 
[1] as a general introduction to Grassmann manifolds and also [5] and its references as an 
application. 

Namely, we treat a hermitian matrix 

( H + \ 

H= G H(n;C) (13) 

\V H ) 

where 

H + G H(k; C), H G H{n- fc;C), V G M(n — fc, k; C) 
for 1 < k < n — 1. In order to make H a direct sum form we prepare a unitary matrix 

/ lfc -tf \ I (l fc + ZtZ)- 1 ^ 

= 17(Z) = 1 ' J J (ID 

V £ In-* / V + 

where Z G M(n — fc, fc; C) is a parameter matrix. U is a map 

U : M(n-k,k;C) — >■ S?7(n) 

and Z is a local coordinate of the Grassmann manifold Gfc(C n ) defined by 

G fc (C n ) = {P G M(n; C) | P 2 = P, P f = P, trP = fc} 
= {UPorf | 17 G U(n)} 
U(n)/U(k) x U(n — k) 

with P given by 

Note that dinicGfc(C n ) = k(n — k) = dimcM(n — k, k; C). The local parametrization of 
Gk{C n ) is more explicitly given by 

p(z) = u(z)p u(zy = | lk ~ z] ) ( lk ) ( lk " zt | (is) 

Z l n -k I \ On-k I V Z l n -k 



where we have used the relation 



-1 / x / x t 

1* \ [ (u + z^z)- 1 \ ( l k -tf 



Z ln-k J \ (In-* + Z&)- 1 \ Z 1 



n—k 



Let us calculate U\,,HUa 



uIjHUm 



■ (16) 



1* & \ I H + Vt \ / i fc 

h + + zty + v^z + zt#_z yt _ h + z^ + zt#_ - ztyzt 

V - ZH+ + H_Z - ZV^Z H_ - ZV^ -VZ^ + ZH+Z^ 
From this we set 

V - ZH+ + H_Z - ZV ] Z = ^ Z\/ f Z + Zi/ + -iJ_Z-y = 0. (17) 
This is just the matrix Riccati equation. Under the condition we obtain the block form 
U]RU= ( (l k + ZtZ)-V*H + (l k + &Z)-V> 

\ (i n - k + zzr r ^H_(i n _ k + ztfyw 

(18) 

where 

H + = H + + Z ] V + V ] Z + Z ] H_Z, H_ = H - ZV ] -VZ ] + ZH + Z ] . (19) 

How to solve the Riccati equation (|T7|) is not known as far as we know. In fact, it is very 
hard, so we must satisfy by finding some approximate solution at the present time. 

Approximation I 

First, by rejecting the quadratic term we have 

ZH + -H_Z = V. (20) 
This solution is well-known to become 

POO 

Z= e tH -Ve- tH+ dt (21) 
Jo 
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under some condition on H_ and H + . See for example [6]. In fact, 

POO 

ZH + -H_Z = / {e tH ~Ve~ tH+ H + - H_e tH -Ve- tH+ }dt 
Jo 



* Ue tH ~Ve- tH +)dt 
dV ' 



o 



under the condition 



-\e tH -Ve' tH+ ]^ 
V 



lim e tH -Ve~ tH+ = 0. (22) 

t— >oo 



Approximation II 

Next, let us consider another approximation. We assume that n = 2m, k = m and V is 
invertible (V G GL(m; C)). Then, by remembering 

ax 2 + 2bx + c = a(z H — ) = — c H 

a a 

we have 

{Z - H^V^jV^Z + (yt)- 1 ^} = V — H_(V ] )- l H + (23) 
from (fTTj) . Here if we can choose Z as 

Z + (V ] Y l H + e GL(m;C) 

then we have a recursive relation 

Z = H^y 1 + {V- H^y 1 H + } z + ^ ]) _ 1H ^ )-\ (24) 



Now by inserting an approximate solution (jzijl into the equation above we obtain the ap- 
proximate solution 



if 



/ e^-Ve^+dt + {V ] Y 1 H + e GL(m; C) 



or 



H + + / V^e tH ~Ve~ tH+ dt G GL{m; C). 

A comment is in order. We don't know at the present time whether our approximate 
solution is convenient enough or not. 

4 Reduction of Riccati Diagonalization 

In this section we give an explicit procedure of Riccati diagonalization. General Hamiltonian 
is 



H 



hi 


V21 


V31 


■ Vn-1,1 


Vnl 


V21 


h 2 


V32 


■ Vn-1,2 


Vn2 


V31 


V32 


h 3 ■ 


■ Un-1,3 


Vn3 


ra-1,1 


Vn-1,2 


Vn-1,3 ■ 


h n -i 


Vn,n— 


Vnl 


Vn2 


Vn3 ■ 


Vn,n—1 


K 



\ 



e H(n;C) 



(26) 



and we write as 

/ h+ yt 

H=\ |, V = {V n i,V n2 , ■ ■ ■ ,U n ,n-l) 

for simplicity. We prepare a unitary matrix 

;i B _! + ztz)- 1 / 2 



(27) 



l„-i 
Z 1 



i + zzt) 



'1/2 



(28) 



where Z = (^i, z 2 , • • • , z n -i). Then the Riccati equation is 
ZV^Z + ZH+ - h n Z - V = <^ 



n—l \ n— 1 

j=l ) 3=1 

for 1 < k < n - 1. (29) 



Note that to solve the equation(s) above explicitly is very hard, so in general we must 
satisfy by constructing some approximate solution. 
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If we can solve the equation(s) then 



U ] HU = 



(l n _! + ZtZ)-V2/ /+ (i n _ 1 + ztZ)-V2 



n— 1 n— 1 n— 1 

= - y](^nj + C.C.) + ^ Z j( H +)jk z k 

3=1 j=l k=l 

and the procedure is reduced to the calculation of 

(l n _! + ZtZ)- 1 /^ + (l n _ 1 + ^Z)- 1 / 2 

, so we must calculate the term (l n _i + Z'l'Z)^ 1 / 2 exactly. 
We write 

Z= {Z 1 ,Z 2 ,--- ,Zn-l) =Z 1 (1,W 2 ,--- ,W n -i) =Zi(l,W), Wj =Zj/ 



ZZ ] = \ Zl \ 2 (l + WW*) = 

3=1 



n-1 

2 



\ Z j\ 



for simplicity. Then 

l n _ 1 + ^Z = l n _ 1 + \z 1 \ 2 
and a unitary matrix given by 



1 W 



gives 



1 W (l + w^w^t)- 1 / 2 

^ ln-2 / I (l n _ 2 + ^t^)-V2 



fir | 1 W^^-r 1 ^ 
o„_ 2 i 1 + ww] [ WW 
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Therefore 



ln-i + tfZ = l n ^ + \ Zl \ 2 (l + WW^U [ 1 ) 



>n-2 



n-1 



3=1 V o 



n-2 



L n-2 



and we have 



-1/2 



(l n _! + ZtZ) 
As a result the reduced Hamiltonian is 



= C/ 



ln-2 



(l n _! + Zt^)- 1 /2 ff+ (l n _ 1 + ^Z)- 1 / 2 
/ 1 \ / 



\ 



ln-2 / 



V 



ln-2 



and we obtain 



V 



Wh+u 



K-2 ) 



( 



\ 



l n-2 



(32) 



(33) 



\ 



(34) 



We have only to continue the reduction process one after another. 

A comment is in order. Note that it is not easy to calculate (l n -k + Z^ Zy 1 / 2 for 
Z e M{k, n - k; C) and 2 < k < n - 2. 
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Note There is no need to calculate (l n _2 + W^W) in U because 



U 



\ 



1 —W 
1 —W 

l n _ 2 



ln-2 / 

/ 



(l+WH/t)^i+ E -i| 2j |2 / 1 W 

(l n _ 2 + W^-i J \ -W/t i n _ 2 

1 \ 



(l+HTTtl^r'^p 



V 



ln.-2 — 



1 w 

l n _ 2 



L n-2 1 

from the definition of U. This is important. 

Last, let us make a comment. If n = 3, namely Z = (zi, z 2 ) then we have a direct method 
(without using W). For 



1 2 + Z^Z 



l + |2i| 2 Z X Z 2 
Z\Z 2 1 + |-22| 



we set 



U 



^1 — ^2 



^2 2l 



G 5C/(2). 



Then we have 



and 



(h + z^zy 1 / 2 = u 



I 



\ 



ki + ^2 



Zi -z 2 
z 2 Z\ 
( M 2 



Vi+kiP+k 2 | s 



v 



+ \ Z 2\ 



Z\Z 2 



Zl z 2 
-z 2 Z\ 

\ 



^1^2 



This form looks smart and will be used in the next section. 



(35) 



/ 
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5 Example 

In this section we apply our method to the following important example 



H 



Since 



#4 



hi a 
a ho 



a h 2 7 
V P 7 h J 



, V = 03, 7 ), Z = (z 1 ,z 2 ) 



the Riccati equation is 



or 



/3zf + 2 2 + (hizi + az 2 ) - h 3 zi — (3 — 0, 
jz% + /3z l z 2 + (a^i + h 2 z 2 ) - h 3 z 2 -7 = 

f3zf + (/ii - /i 3 )2i - (3 = -z 2 (a + 721), 
74 + (/i 2 - h)z 2 - 7 = -zi(a + /fofc)- 



Let us solve the equations. We get 24 by solving the equation 
{Pl(h - h 2 ) - a(3 2 + a^}z\ + 

M(h - h2)(h - h 3 ) + 2\a\ 2 - \(3\ 2 - | 7 | 2 } - a(3(h + h 2 - 2h 3 )]z 2 + 

- h 3 )(h 2 - h 3 ) - \a\ 2 - \I3\ 2 + 2| 7 | 2 } + /3 7 (-2/n + /i 2 + /i 3 )]*i + 
/3 2 7 + a/3(/i 2 - /13) - « 2 7 = 

by use of Cardano formula and z 2 by solving the equation 

7z 2 + (h 2 — h 3 + fizi)z 2 + az\ -7 = 0. 
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If (zi,z 2 ) is given then the reduced Hamiltonian becomes 

\-l/2 



(h + tfzy^H^h + tfzy 
( 



1 



2^ 2 



+ 1^1 



Z\Z2 



Vi+kiP+i^i 2 

^ 



— / — Zi Zo 
^ +kl| 2 



\ 



x 



a 



+ jzi + /3z 2 + hsZ\Z2 
a + jzi + f3z 2 + h^zxz 2 h 2 + "yz 2 + ^z 2 + h 3 \z 2 \ 2 



kil s 



+ *2 



■Z1Z2 



Zl£2 



y/l+\zi\ 2 + \z2\- 
_^_^ = _ ^ ^ M 2 

V Vi+I^l'+I^l 2 1 2 y/i+\*i\ 2 +M' 

Therefore we have only to solve the equation 

(z 2 + (An - k 2 )z - C = 

as shown in the example in section 2. 
Finally we obtain three eigenvalues 



- z x z 2 



+ l*i I 



\ 



Ai 
A 3 



fci + (z + + k 2 \z\ 



^2 



&2 - C z - C z + h\z\ 



h - (13 zi + /3zi) - (-yz 2 + ^z 2 ) + hi\z x \ 


2 + aziz 2 + otZ\Z% + /i 2 *2 2 


1 + 


*i 


2 + 


1*2 


2 



(40) 



(41) 



(42) 



under (|38l), (|39|) and (j4T|). 

In the process of calculation MATHEMATICA or MAPLE is indispensable (calculation 
by force is very hard). 



6 Discussion 

In this paper we presented a geometric approach to diagonalization of hermitian matrices. 
The advantage of our method is quick diagonalization, while to obtain eigenvalues is left in 
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the final step. Our method is in a certain sense reverse process of the standard one, so they 
are dual each other. 

It is not clear at the present time whether our method is convenient enough or not. 
Further work will be needed. 

Last, let us make a comment on "Riccati structure" of Quantum Mechanics, see for 
example [Zj. We consider the harmonic oscillator given by 



In order to solve the model we usually define the annihilation and creation operators like 





Then it is well-known that 




a) a = H 



1 



aa 



2 



Next we define another annihilation and creation operators like 




with unknown (3(x). 

Then, in order to satisfy the relation btf 



H + | /3(x) must satisfy the equation 



(3' + (3 2 



1 + x 2 . 



This is a special type of Riccati differential equation. 



Appendix 



In the appendix we consider the important example 



( h 



i 



a 




H 



a h 2 



7 



V P 7 



^3 / 
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once more. We want to diagonalize the matrix at a time. 

For the purpose we prepare a unitary matrix coming from the flag manifold of the second 
type (in our terminology, see [8], [9] and also [ID], p]) SU(3)/U(1) x £7(1) 



U = U(x,y,z) 



1 



-(x + yz) xz — y 



x Ai — x(x + yz) —z 
V V zAi - y(x + yz) 1 j 



( 



v / A7A^ 



= U M U D 



where Ai and A 2 are given by 



\ 



(43) 



Ai = 1 + \x\ 2 + \y\ 2 , Ai = 1 + \z\ 2 + \xz - y\ 2 . 
Let us calculate U^HU. Since U^HU = UdU\ { HUmUd we have only to calculate 



W m HUm because Up is diagonal. Namely, 



UlHUM 



( 1 x y 

— (x + yz) Ai — x(x + yz) zA 1 — y(x + yz) 
y xz — y —z 1 

( 1 — (x + yz) xz — y "\ 

x Ai — x(x + yz) — z 

V y zA x - y(x + p) 1 y 

/ IDu W 2 1 U>31 ^ 

u>21 W22 W32 

V W 3 i W 32 ^33 / 



\ / hi a P \ 

a h 2 7 

/ V p 7 ^ / 



(44) 
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where 

W21 = —(x + yz)(hi + xa + y/3) + {A 1 - x(x + yz)}(a + xh 2 + 2/7) + 

{2A1 - y(x + yz)}(/3 + ary + yh 3 ), 
w 3 i = (xz - y)h\ - za + (3 + {(xz - y)a - zh 2 + 7)2; + {(xz - y)/3 - 27 + /i 3 }y 

= (xz - +0:0; + y/3) - z(a + x/i 2 + 2/7) + (/? + ary + y/i 3 ), 
W32 = — + yz)[{(xz - 2/)^i - za + /?} + {(xz - y)a - zh 2 + 7}x + 

{(xz - y)j3 - 27 + h 3 }y] + Ai[(xz - y)a - zh 2 + 7} + {(xz - y)/3 - zj + h 3 
= -(x + yz)w 3 i + Ai[(xz - y)a - zh 2 + 7} + {(xz - y)(3 - 27 + fo 3 }z]. 

Here by setting u> 2 i = w 3 i = w 32 = 0, we have 

= -(x + yz)(hi + xa + y/3) + {Ai - x(x + yz)}(a + xh 2 + y^y) + 

{zAi - y(x + yz)}((3 + 2:7 + y/i 3 ), 
= (xz — y)(hx + xa + y/3) — z(a + xh 2 + y*f) + (f3 + x^ + yh 3 ), 
= (rrz — y)a — zh 2 + 7 + {(xz — y)(3 — z^ + h 3 }z. 

Some calculation gives 



x(h\ + xa + y[3) — (1 + 




2 )(a + x/i 2 + 2/7) + + ^7 + |//i 3 ) 


+ + + xy(a + rc/i 2 + 2/7) - (1 + 




2 )(/3 + :r7 + 2//i 3 ) 



2/(/ii + xa + - (/3 + X7 + y/i 3 ) 

.2- — — 

x(h\ + i« + y[3) — (a + x/i 2 + t/7) ' 
(x/3 — 7)-2 2 + (^3 — ^2 + xa — y(3)z — ya + 7 = 

in terms of Ai = 1 + |x| 2 + \y\ 2 . 

It is not easy at the present time to solve the equations at a time. 
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